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ABSTRACT 

An all digital simulation of the existing boiler and control system 
of CVA-63 was carried out. A significant savings in computer time over 
previous simulations was realized. 

The parameter-plane method was used to search for new operating 
points for each of the system minor control loops. Individual loop 
transient response was tested with a simulation program to demonstrate 
improved response with new operating points. 

The entire system response was then investigated to demonstrate the 
overall effects. Recommendations were made for new controller settings 
along with component variations in order to improve over-all system 
response . 
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I. INTRODUCTION 



Naval boilers differ from most other stationary type boilers in 
that they are called upon to answer widely varying load demands in a 
minimum amount of time. Feedwater, fuel oil and combustion air must be 
supplied in varying, but accurate, amounts to maintain sufficient steam 
flow as required by ship operation. There are maneuvering conditions 
which might require boiler load to change from ten or twenty percent to 
over one hundred percent in a matter of a few seconds. 

Automatic combustion control systems have been developed by various 
manufacturers to control the three inputs; oil, air and water. Most of 
the ACC systems up until the present have been of the pneumatic variety; 
i.e., pneumatic sensing and valve control. Individual components of 
systems of this type range from simple bellows type sensors to the ex- 
tremely complex masses of springs, levers, valves, gages and nozzles 
found in pressure relays and transmitters. 

This thesis will investigate the automatic combustion control system 
including a main boiler as found on USS KITTY HAWK (CVA-63) . An all- 
digital computer model will be developed for testing and analysis purposes. 
Each of the three major control loops, fuel flow, air flow and feedwater 
flow, will be investigated independently and the effects of any changes 
or improvements will be tested in the system simulation. No radical 
changes to the overall system will be proposed but readjustment of vari- 
able parameters within the existing system will be stressed with emphasis 
on improved overall system performance. 

A # THE BOILER 

The boiler to be modeled is a Foster-Wheeler D-type as found on 
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CVA-63. The boiler consists of an upper steam drum connected by tubes to 
the lower water drum, screen wall and waterwall headers, A vertical super- 
heater is integral to the boiler. This unit is capable of generating 
318,000 lb. /hr. of steam at 120 percent overload conditions. Cruising 
conditions are considered to be 117,500 lb. /hr. of steam which is 44.7 
percent of full boiler load. See Appendix A and Ref. 1 for more design 
and performance data. Two boilers are operated together and controlled 
by one ACC unit due to physical location on shipboard. Boiler transfer 
functions were obtained from Ref. 2 and were valid for cruising condi- 
tions only. For this reason, small variations in boiler load were used 
and the model developed was of the perturbation variety. When data on 
other conditions of load become available, the model may be easily modi- 
fied to suit the design engineer. 

B. AUTOMATIC COMBUSTION CONTROL 

The automatic combustion control system encompasses two classifica- 
tions of pneumatic control, combustion control and feedwater control. 

While the feedwater portion is not essential for combustion, it is usually 
included in any discussion of ACC systems since it is very necessary for 
operation of the overall system. 

The primary function of the combustion control portion of the system 
is to maintain constant boiler superheater outlet pressure under all load 
conditions. It must also insure that fuel and air are mixed in the 
proper proportions to maintain optimum combustion efficiency. Normally, 
a safety feature is p nvided with the system whereby, in the event of loss 
of control air pressure, all valves return to a neutral position which 
reduces boiler firing to a stable minimum. 
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FIGURE 1(a) 
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FIGURE I (c.) 
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FEEDWATER FLOW CONTROL 




FIGURE 1(d) 
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CONTROL SYSTEM 
INTERCONNECTION DIAGRAM 

FIGURE 2 
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The feedwater control portion of the system must maintain a balance 
between feedwater flow to the boiler and steam flow from the superheater. 
It must also maintain the water level in the steam drum at a normal 
steaming level under all conditions of boiler load. This level control 
is a most important function since high water level can mean possible 
liquid-state carryover into rotating turbines and subsequent destruction 
of blading while low water can expose the ends of tubes within the boiler 
leading to sagging and weakening due to extreme heat without adequate 
circulation cooling* Level control is complicated by the phenomena of 
shrink and swell. When boiler load is increased rapidly, large amounts 
of steam are removed from the boiler and the water in the steam drum will 
surge to a high level. The reverse is true when the load is suddenly 
decreased. The thermodynamic explanation of shrink and swell is to be 
found in Ref. 3 while more discussion is found in this paper in Section 
II. A. 6. 



C. OPERATION OF AUTOMATIC COMBUSTION CONTROL SYSTEM 

The steam pressure transmitter (Fig. Id) must transmit a linear 
pneumatic signal directly proportional to pressure variations in the two- 
boiler main steam header. At the normal operating pressure, 1200 psig., 
the output of the steam pressure transmitter is 9 psig. The steam pres- 
sure controller is of the set-reset variety, also known as proportional 
plus integral control. When the steam pressure signal from the steam 
pressure transmitter decreases, the output of the steam pressure controller 
increases proportionally and conversely. The proportional band setting 
of the controller determines magnitude of output variation for a given in- 
put and the reset needle valve of the controller determines the rate 
response. 
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The output of the steam pressure controller goes both to the air 
flow controller and to the low-pressure signal selector. This output 
becomes the demand signal for the air flow control loop. The air flow 
controller compares this signal to the actual air flow signal which is 
proportional to the flow of combustion air across the boiler air regis- 
ters. The air flow controller is also of the proportional plus integral 
type and its output signal is sent to the airflow ratio relay and damper 
actuators. There are two factors which affect the quantity of air flow, 
damper position and forced draft blower turbine speed. The ratio relay 
preceeding each blower controller insures parallel operation of blowers 
for there can be one, two or three blowers in operation for each boiler 
at any one time. It has been found that the use of dampers is not only 
unnecessary but that they have actually contributed to some air flow loop 
instability problems. Dampers were completely eliminated on USS AMERICA 
(CVA-66) (See Ref. 4). 

The air flow transmitter senses pressure differential between wind- 
box and furnace and transmits a signal proportional to this difference. 
This signal, representing actual airflow across the burner registers, is 
sent to a ratio relay where the important air-fuel ratio may be adjusted 
for optimum combustion efficiency. 

The output of the ratio relay goes to the air flow controller already 
discussed and to the low-pressure signal selector. This signal selector 
transmits the lower of its two input signals, boiler load demand from 
the steam pressure controller and air flow signal from the fuel-air ratio 
relay. The lower of these two signals is used as a demand signal to the 
oil flow portion of the system so that, under increasing load, oil sup- 
plied to the furnace can never exceed the required air flow for proper 
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combustion. This prevents the formation of black smoke in the furnace 
when the load is suddenly increased. 

The output of the signal selector goes to a characterizing relay 
which converts it to a corresponding non-linear signal conforming to the 
burner curve of return oil pressure vs, burner oil rate. Fuel oil is 
supplied at a constant supply pressure and the amount burned in the 
furnace is controlled by closing down on the return valve in the recircu- 
lating system. 

The feedwater portion of the control system consists of a thermo- 
static drum level sensor, a combining relay and the feedwater regulating 
valve . 

The level sensor consists of a stainless steel expansion tube in- 
clined across the steam drum with its center at normal steaming water 
level. The upper tube end is rigidly fixed and the lower end tube move- 
ment is amplified through a lever system attached to a pilot valve stem. 
Level sensing motion comes from contraction due to temperature difference 
between saturated steam temperature (upper half) and boiler water tempera- 
ture (lower half). The pilot valve converts the mechanical motion into a 
proportional air output signal. 

A boiler steam flow signal and water flow signal are combined with 
this level signal in the combining relay using a combination of diaphragms 
and levers. The resultant signal is a total water demand signal and 
controls the feedwater regulating valve. The feedwater valve positioner 
controls the flow of feedwater by converting the demand signal to mechani- 
cal valve motion. 

Now that a general system description has been given, consider the 
effects of changing the boiler load. 
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1 * System Action During Increasing Boiler Load 

The first noticable effect on increasing boiler load is a drop 
in superheater outlet pressure. This is reflected by a decrease in the 
output of the steam pressure transmitter which causes a corresponding in- 
crease in the output of the steam pressure controller. 

This represents an increasing air demand signal which is greater 
than the actual airflow signal. This causes an increase in blower speed 
and opens dampers further until flow matches demand. 

When load increases, the signal transmitted by the low-pressure 
signal selector is the smaller air flow signal* As the air flow return 
system receives an increasing signal, the fuel oil return valve is closed 
proportionately causing more oil to be burned in the furnace. This, in 
turn, causes an increase in steam pressure until equilibrium is once again 
reached. 

The effects of boiler water swell are reduced by the action of 
the combining relay. As load increases, the steam flow signal increases 
to a new, higher fixed position. The swelling water level produces a 
false signal (actually, more water is needed, not less as indicated at 
first) which is balanced by the steam flow signal. Relay settings deter- 
mine sensitivity. As the level swell boils away, the level signal will 
decrease to its true value and a demand will be made for more water. As 
water flow increases to a higher value, its signal will offset the steam 
flow signal and equilibrium will be reached once again. 

2 . System Action During Decreasing Boiler Load 

When boiler load is decreased, the superheater outlet pressure 
increases causing an increase in steam pressure signal from the steam 
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pressure transmitter. This increasing signal causes a proportional de- 
crease in steam pressure controller output. Since air flow is now greater 
than necessary, the output of the minimum signal selector becomes the lower 
steam pressure signal, decreasing fuel oil demand. The air flow response 
lags behind at this point, providing an excess of combustion air until 
blowers slow and dampers close. 

The decrease in fuel oil demand opens the return valve and the 
pressure returns to equilibrium due to the smaller quantity of fuel oil 
being burned in the furnace. 

The water flow portion will act oppositely to the manner in 
which it acted for a load increase. Decreased steam flow demand will be 
partially (or fully, depending on adjustment of combining relay) offset 
by shrink in boiler water level signal. As this shrink effect dies out 
and the drum level starts to increase, the output of the combining relay 
will increase and the feedwater regulating valve will close. As water 
flow decreases, the water flow signal offsets the steam flow signal in 
the combining relay and equilibrium is reached. 

D. IMPROVEMENTS NEEDED 

Control of a boiler and its associated auxiliary equipment is compli- 
cated considerably by the interaction present among its control variables. 
The representation of the boiler and control system in Figs. 1(a) through 
1(d) and the interconnection diagram of Fig. 2 demonstrate the coupling 
which exists. All three major system variables; steam flow, oil flow, 
and water flow affect both boiler steam pressure and boiler water level 
in varying ways. The transfer functions seem complicated enough at first 
inspection but when one stops to consider that these are valid for cruising 



22 



conditions only and that at varying boiler loads they become dynamic 
functions of the load and rate of load change the control problem becomes 
ove rwhelming. 

Ideally, if the interaction between inputs and outputs could be re- 
moved or balanced in some way, the control problem would be greatly 
reduced in scope. In the case of the boiler, physically forcing non- 
interaction is impossible. The alternative then, is to design a non- 
interacting controller. This has, in fact, been done on paper (See Ref. 

5), but the controller transfer functions which result for this system 
are not physically realizable and could not be simulated using the digital 
simulation language to be discussed in the next section. See Appendix B 
for design information on the non- interacting controller. 

Any improvement then, in control of this boiler, is limited to 
complete redesign of the entire system, or optimization of the existing 
system. 

Alternative control systems have been proposed. Digital control 
(Refs. 6 and 7) and analog fluidic control (Ref. 8) are two of the better 
known alternatives. The first was investigated by Naval Electronics Labor- 
atory Center but this investigation was terminated in June 1968 due to 
unfavorable reaction by Naval Ship Engineering Center*. Analog fluidic 
control is under test and development presently at NAVSEC, Philadelphia. 

It is safe to assume that this testing and design will come to fruition 
in a new, fast, state-of-the-art control system for ships of the future. 

What are we to do with the ships of today? Economics forbids the 
complete replacement of several hundreds of ACC systems on ships which 
have a life expectancy of between twenty and thirty years. Hence, let us 
look into the possibility of improving present systems at a minimum of 
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expense. At the very least, we can hope to explore methods which will 
enable the future ship control designer to design better systems. 

Exactly what is needed in the way of improvements? Basically, it 
is a faster responding system with less overshoot of variable response. 

For example. Ref. 9 states that specifications on several ships required 
boiler water level to remain within + 4 inches of set point during linear 
changes in boiler load from ten percent to eighty percent in 45 seconds. 
Subsequent shipboard operation produced high water conditions out of 
sight in the boiler gage glass. This was due to maneuvering the ship with- 
in the load specifications but in only 10 seconds instead of the specified 
45 seconds. Although the system would probably meet the original specifi- 
cations, response could probably be improved somewhat through a judicious 
"tuning 11 of the control system. 

This thesis investigates the CVA-63 system and makes recommendations 
for readjustment of controller settings or replacement of component parts 
as necessary for better (faster) overall system response. The parameter 
plane method of optimization is used and it is hoped that even if system 
response is not dramatically improved, the methods used will prove en- 
lightening for others who may follow. 
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II. THE SIMULATION 



The boiler and automatic combustion control system were simulated 
utilizing the IBM developed Digital Simulation Language (DSL). DSL is 
a non-procedural, problem-oriented language which will accept problems 
in ordinary differential equation form or as analog block diagrams. Since 
an excellent block diagram of the CVA-63 system was available (Ref. 2), 
this program provided a relatively simple, efficient means of simulation. 
References 10 and 11 describe the DSL program in more detail. The vari- 
ables indicated in parenthesis refer to the system diagrams, Figs. 1(a) - 
1(d) and the interconnection diagram, Fig. 2. 

A. BOILER 

The seven transfer functions appearing in Fig. 1(a) are the result 
of tests conducted at NBTL (Ref. 2). These functions are the result of 
a combination of both pulse testing (Ref. 12) and of periodic input fre- 
quency response analysis. 

1. Steam Pressure / Steam Flow (PDS/GS ) 

This response essentially represents a mass balance situation 
and is accurately represented by a straight integration at cruising condi- 
tions . 

2 . Steam Pressure / Fuel Oil Flow (PDF/GF ) 

This can be represented by a first-order lag, first-order lead, 
and a straight integration, all series connected. 

3. Steam Pressure / Water Flow (POW/GW ) 

A second-order lag (complex poles) and series connected integra- 
tion represent the effects of feedwater flow on steam pressure at cruising 
conditions . 
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4. Water Level / Water Flow (LW/GW ) 



As in 1. on the preceeding page, a straight integration ac- 
curately represents this simple mass balance situation. 

5 . Water Level / Oil Flow (LF/GF ) 

Pulse tests by NBTL yielded a transfer function consisting of 
integration and series-connected lag-lead. 

6. Water Level / Steam Flow (LST/GS ) 

The response of water level to steam flow rate is shown by the 
parallel combination with outputs LSA and LSB in Fig. 1(a). Water level 
in the steam drum should integrate proportionally to the mass balance of 
steam flow-water flow difference. This integration, however, is of low 
natural frequency and is initially overshadowed by the shrink/swell 
phenomenon discussed previously. Fig. 3 gives the response, LST, to a 
step input. Fig. 4 give'' the same response to the triangular pulse input 
used for response tests in the system (Fig. 5). It can be seen that the 
effects of the integration are not present for the first several seconds. 

B. FUEL OIL FLOW CONTROL 

Also called steam pressure control loop, the fuel oil flow control 
loop is a closed control loop consisting of the following elements; 

1 . Steam Pressure Controller (PZP/PH ) 

This is a typical proportional plus integral controller. Values 
for proportional and integral gain were measured by NBTL using sinusoidal 
testing techniques and are shown in Fig. 1(b). These gains are adjust- 
able and will be utilized to set the loop to secure "optimum" performance. 

2 . Low Pressure Signal Selector (PM, PQR/PSF ) 

This acts as a switch only and, as such, has no measurable trans- 
fer function. It is conveniently simulated in DSL by using a Function 
Switch (see computer output). 
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FIGURE 3 
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WATER LEVEL RESPONSE TO A STEP INPUT 




FIGURE 4 
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WATER ... LEVEL RESPONSE TO A PULSE INPUT 



3. 



Characterizing Relay (PSF/PRF ) 



This has a two-pole response as determined by NBTL using the 
periodic input method. 

4. Fuel Oil Control Valve (PRF/PF ) 

A quadratic lag function was obtained for this valve. A small 
time delay (0.023 seconds) is also present due to valve-stem stiction. 

This delay appears in all of the pneumatic/mechanical valves and is present 
due to the physical limitations of the valves themselves, and their in- 
ability to respond instantaneously to a demand. 

5. Steam Pressure Transmitter (POUT/PP ) 

Analysis of transient response tests by NBTL yielded a quadratic 
lag transfer function as shown in Fig. 1(b). 

C. WATER LEVEL/FLOW CONTROL 

The water control loop is a M three-element M system, responding to 
changes in water level, water flow, and steam flow. Individual elements 
(see Fig. 1(c)) are as follows: 

1. Feedwater Control Valve (PDW/GW ) 

The transfer function for this valve was found by NBTL to be 
third order, consisting of a first-order lag, attributed to the capaci- 
tance of the pneumatic diaphragm chamber, and a second-order lag due to 
the valve positioner. A 0.45-second pure time delay was thought to be 
due to valve-stem friction in the packing gland. 

2. Combining Relay (PSA ,PWA/PSW ) 

Dynamic measurement could not be made by NBTL and this com- 
ponent is approximated by pure gains. 

3 . Thermostatic Drum Level Indicator (L/PL ) 

At cruising condition loads, the thermostat may be represented 
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by a pure gain. NBTL was unable to make dynamic measurements due to 
limitations in the generation of a suitable water level input signal, 

D. AIR FLOW CONTROL 

The air flow loop is a typical feedback control system and the 
individual component transfer functions were derived by NBTL and are as 
follows : 

1. Controller (POA/PEQ ) 

This controller is a typical set-reset controller with adjust- 
able integral and proportional gains. 

2. Ratio Relay (PG/POA and PQR/PQA ) 

The transfer function obtained by NBTL is as shown in Fig. 1(d). 

3 . Forced Draft Blower Actuator and Valve (GB/PQ ) 

The actuator and valve exhibit characteristic underdamped 
response and time delay which is normal with components of this type. 
Valve-stem friction is not as pronounced as in the feedwater control valve 
due to decreased valve-stem travel and the fact that the medium being 
controlled is steam instead of water. 

4 . Main Forced Draft Blowers (QA/GB ) 

Direct frequency response testing was used by NBTL to derive 
the transfer function shown in Fig. 1(d). 

5. Air Flow Transmitter (PQA/QA ) 

This component is represented by a quadratic transfer function 
and appears in the feedback system of the air flow control loop. 

E. INITIAL SYSTEM TESTS 

Individual components were simulated and interconnected in DSL (see 
program listing) to produce a system simulation for transient response 
testing. Since transient response data for this boiler and combustion 
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control system are almost non-existent, the author had to rely on personal 
experience as Boiler Officer on USS INDEPENDENCE (CVA-62) from 1964 to 
1966, General shape and accuracy of system response were also verified 
by J. W. Banham, Jr., head of the Machinery Automation Systems Depart- 
ment of the Naval Ship Engineering Center, Philadelphia .Division. 

Inspection of the completed simulation program on page 96 will re- 
veal the manner in which individual components and loops are interconnected. 

It was found necessary to put limitations on certain pneumatic signals 
to simulate the loading effects of a controller or relay being overdriven. 
This portion of the simulation appears at the end. 

A triangular shaped input was selected for transient response testing. 
Since the model is of the perturbation variety, the use of a step input 
was deemed impractical and unrealistic. The triangular shaped input seemed 
to be the closest meaningful input signal which most nearly resembled a 
true impulse. The input is shown in Fig. 5. 

System responses to this input are shown in Fig. 6 through Fig. 10. 

It must be remembered that variations are measured from a steady cruising 
condition value; i.e., a steam pressure variation of -8 psig is actually 
1192 psig drum pressure, 8 psig below the nominal 1200 psig. 

The value of a component by component simulation is readily apparent 
at this stage. Any intermediate signal may be designated as output for 
analysis purposes. Components which are overloaded or driven into satura- 
tion are easily discovered and individual components may be replaced with 
new or improved components to test the effects on the overall system 
response . 

The effect of the underdamped feedwater control valve on water flow 
rate is readily apparent in Fig. 6. Later transient response tests show 
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the effects of varying time delay in the transfer function and of 
changing one of the gains in the combining relay. 

The complete system simulation shows a great time saving over pre- 
vious all-digital simulations (Refs, 7 and 13). Run time for the DSL 
program has a 2:1 correspondence with real or problem time. For example, 
two minutes of simulated boiler operation takes one minute of execution 
time on the IBM 360/65. While this depended somewhat on the size of the 
integration step and on the number of variables output, previous simula- 
tions were approximately ten times slower by comparison. 

Simulation of individual loops was used for transient response testing 
but the simple programs used are subsets of the overall system program 
and do not merit further discussion. 
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Ill, INDIVIDUAL LOOP ANALYSIS 



Since system complexity precluded mass analysis of the entire ACC 
system at one time, it was decided to separate the system into three 
minor loops. Obvious choices for analysis were the three major system 
variables, combustion air, fuel oil and feedwater. Although interaction 
is present throughout this system, it was ignored for individual loop 
analysis in the hope that improvements made in each loop would be signifi- 
cant in the whole system. This obviated the massive problem of trying to 
counter the effects of interaction to yield accurate loop responses. 

Several analysis methods exist for determination of controller para- 
meters to permit "optimum" performance. Reference 9 gives an excellent 
overview of individual component relationships to system performance. 
Essentially however, due to the demands of shipboard maneuvering, the 
performance requirement can be easily stated: 

The ACC system must be capable of maintaining stable operating 
conditions under widely varying load conditions with a given 
set of machinery (boiler, blowers, valves, etc.) and it must 
make transitions between operating points in a minimum amount 
of time. 

This requirement is more easily stated than met. The designer is limited 
to passive components and more precisely, to the relative slow pneumatic- 
mechanical components which are presently available. 

Root-Locus, Bode, and Black-Nichols are among analysis methods 
presently used in the design and analysis of systems of this type. 
Reference 14 is the first known to apply the parameter plane technique 
to the minor control loop of a boiler ACC system. 
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The parameter plane is a two-parameter method which determines the 
values of two parameters (0Cand(3) that are required to insure that the 
characteristic polynomial has roots at a prespecified location in the s- 
plane. By choosing a suitable line on the s-plane and mapping it through 
certain algebraic transformations (see Ref. 15) one may obtain a curve on 
the Oi “P plane which may be used to select "optimum 11 values of 0(and |3 
to give a desired response. In a dynamic system of this type where speed 
of response was important, three mapping transformations were deemed neces- 
sary. Constant- C0 n lines (U) n is the natural frequency of an approximate 
second-order system) map from circles of constant radius on the s-plane. 
Constant- £; lines ( ^ is the damping ratio of the approximate second-order 
system) map from radial lines on the s-plane. Constant- CT lines map from 
real root locations on the real axis of the s-plane. 

One difficulty encountered was the mathematical representation of pure 
time delay as present in all three minor loops. It was necessary to have 
this delay in a mathematical form to determine the characteristic poly- 
nomial for each of the loops. The Pade Approximation of transport lag is 
among the most widely known and was selected (Ref. 16). Varying degrees 
of approximations were available, ranging from first to twelfth order with 
accompaning increases in accuracy with higher order. Since the use of the 
parameter plane was dependent on the determination of the minor-loop 
characteristic polynomial, a fourth-order approximation was chosen for all 
transport lag representations. This approximation gave a satisfactory 
mathematical result and increased the highest power of the characteristic 
polynomial by four. The general form of the approximation appears in Ap- 
pendix C along with the actual equation used in each of the loops. 
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Since hand calculation of mathematical mapping transformations 
would be prohibitive, a computer program, developed at the Naval Post- 
graduate School under the direction of Dr. George J. Thaler, was used 
for parameter plane graphical results. 

Algebraic manipulation of transfer functions was also performed on 
the computer. FORMAC (FOrmula MAnipulation Compiler), an algebraic manipu- 
lation language (see Refs. 17 and 18) was extremely useful for develop- 
ment of loop characteristic polynomials and should prove useful to the 
control engineer faced with tedious algebraic manipulation. 

All graphs presented in this thesis are the result of computer runs 
and were drawn directly by a CalComp plotter at the computer facility at 
the Naval Postgraduate School. 

A. AIR FLOW LOOP ANALYSIS 

The air flow control loop is a fairly typical feedback system and is 
shown in Fig. 11 which is a simplified version of Fig. 1(d). The re- 
sulting loop characteristic polynomial for the closed loop system is of 
fifteenth order and is shown in Appendix C. 

The parameter plane approach was used and the selected variable 
parameters were the proportional and integral gains of the air flow 
controller. These variables were thought to be the most easily acces- 
sible and most easily adjusted. The proportional gain in the original 
system had been set at 1.67 and the integral gain at 0.024. For para- 



tional gain alpha (CX). The original system operating point is marked 
with an X on Fig. 12, the air flow loop parameter plane. 



meter plane purposes, integral gain was designated 




propor- 
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AIRFLOW LOOP PARAMETER 

PLANE 



FIGURE 12 



In order to demonstrate the accuracy and validity of the parameter 
plane, several runs of transient response were made using a simple DSL 
simulation of the system in Fig. 11. 

It was decided to hold one parameter constant while varying the 
other in order to dramatize the effect of moving in a straight line on 
the parameter plane. In Fig. 13, alpha was held constant at 1.67 while 
beta was allowed to take on values of 0.01, 0.024, 0.1, 0.2. Note that 
the second curve; 0( = 1.67, = 0.024 is the transient response of the 

original system. Input for all transient response tests in this section 
was a unit step. Since motion along a line of constant alpha constituted 
motion along a line of nearly constant zeta (damping), no increase in 
response time was anticipated. Figure 13 bears out this expectation. 

Motion along a line of constant beta however, was nearly perpendicular 
to the lines of constant zeta, and rapid changes in transient response 
were expected. Figure 14 depicts such transient response; holding beta 
constant at 0.024, alpha was allowed to take on values of 1.0, 1.67, 2.0, 
2.5, 3.0, and 3.5. Again, for comparison, the 01= 1.67, ^ - 0.024 

curve is the response of the original system. At 0(= 1.0, which is nearly 
on the £> = 1.0 curve, the system exhibited the characteristic response of 
a critically damped second-order system as predicted on the parameter 
plane . 

While no one point could be selected as the "optimum 11 operating point, 
several qualitative statements could be made. It is obvious from Fig. 13 
that beta, at least within the range of investigation, had little effect 
on the response time of the system. On the parameter plane, if alpha is 
allowed to remain constant, increasing beta will move the operating point 
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along a nearly constant-^ path until beta exceeds about 0.38 at which 
time ^ begins to decrease. The effect of variation of alpha was more 



readily apparent. Rise times were considerably improved with alpha values 
of 2.0 and 2.5 while overshoot was not considered to be excessive (Fig. 
14). On the basis of this qualitative analysis, the ’'improved 11 operating 



a G on the parameter plane. Fig. 12. Section IV demonstrates the response 
of the entire system with this new operating point. 

B. OIL FLOW LOOP ANALYSIS 

Since the position of the low-pressure signal selector determines 
the configuration of the oil flow control loop, a choice had to be made. 
During an increasing load, the air flow system lags behind and it is the 
output of this loop (PQR in Fig. 2) which becomes demand for oil flow. 
Hence, the air flow loop is an integral part of the oil flow loop when 
load is increasing. Increasing load conditions, before the air flow loop 
responds, the demand signal (PM) from the steam pressure controller is 
lower and therefore becomes the demand signal. 

It was decided to try to improve the loop in the simpler of the two 
configurations since combination would have led to a characteristic poly- 
nomial of the twenty-eighth order. It was hoped that improvements made 
to the air flow loop would help the response when this loop was a part 
of the oil loop during a load increase. 

The simplified oil loop diagram appears in Fig. 15. The resulting 
characteristic polynomial of thirteenth order was determined using the 
FORMAC program and is found in Appendix C. 



point of 0( = 2.5, 




This point is marked with 
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FIGURE 15 



SIMPLIFIED OIL FLOW 
CONTROL LOOP DIAGRAM 



As was done with the air flow loop, the parameter plane was plotted. 
This time, the variables of the steam pressure controller were used for 



parameters. Alpha was chosen to be the proportional gain and beta the 
integral gain. The resulting graph is shown in Fig. 16. 

Several points on the parameter plane were chosen for DSL simulation 
testing. With alpha held constant at 6.67, beta was allowed to take on 
the values 0.001, 0.014, 0.1, and 1.0. These curves are shown in Fig. 17. 



system as marked on the parameter plane with an X. Three curves are very 
close together, the result of moving vertically along a nearly constant- 
line as before. The fourth curve, p - 1.0, represents a point off the 
parameter plane but close inspection will reveal that the value of 
has been drastically reduced for this operating point. The transient 
response bears out this hypothesis, demonstrated by the large overshoot 
of a characteristically underdamped system. Again, with beta held con- 
stant at 0.014, alpha was allowed to vary from 4.0 to 10.0. These re- 
sponse curves are shown in Fig. 18. The critical damping ratio, £> = 1.0, 
was not approached with the oil loop but extrapolation of the curves of 
Fig. 16 to an alpha of 2.5 will lead to the expected response. The 
original system response is shown for comparison as before ( # = 6.67, 



Again, it was difficult to make a quantitative analysis to determine 
an "optimum" operating point. It must be remembered that the results 
were valid only for an increasing load, the loop being considerably 
different during a decreasing load. A new operating point was selected 



Note that 0( = 6.67, 




point of the original 



P » 0.014). 



on a qualitative basis at d = 8.0, 
with a 0 on the parameter plane. 




point is marked 
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C. WATER FLOW LOOP ANALYSIS 

The water flow control loop is shown in simplified form in Fig. 19* 

In order to maintain true closed-loop action for the resulting analysis, 
the loop had to be reduced as shown from the loop in Fig. 1(c). 

The integration in the feedback circuit made the system character- 
istically type one. Hence, response to a step input contained no steady- 
state error, and for unity input, unity output was obtained. As was found 
in the original system response (see Fig. 6) the response of the water 
loop left something to be desired as far as oscillatory action was con- 
cerned. 

The loop characteristic polynomial was determined using the FORMAC 
program and is given in Appendix C. 

Choices for parameter plane variables were somewhat limited in this 
case. Since there was not a proportional plus integral controller present 
as in the oil and air loops, other variables had to be selected. A 
simple inspection of the loop of Fig. 19 revealed two likely parameters. 

The gain of the thermostat (2.37) was selected as alpha and the gain of 
one of the combining relay factors (0.014) was selected as beta. 

Effects of the rather large (0.45 sec) time delay were questioned 
in Ref. 2, so several transient response runs were made of the entire 
system and results are shown and discussed in the next section, IV. D. 

The parameter plane plot is shown in Fig. 20. Note that the operating 
point of the original system, # = 2.37, |3 « 0.147, which is marked with 
an X, lies outside of the ^=1.0 curve. This indicated that the dominant 
roots were not complex, but real. A constant-sigma, or real root, line, 

<r = 0.52, passes through the operating point. This indicated that one 
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of the dominant roots lies at <T » *0.52 on the root locus plot. In order 
to verify this result, a computerized root-locus plotting program, also 
developed at the Naval Postgraduate School, was used. 

Since the root locus is a one-parameter method, it was decided to 
fix alpha (thermostat gain) at 2.37 and more vertically on the parameter 
plane along a straight line and vary beta from 0.1 to 1.0. Since =* 1.0 
appears to be above the £> = 0.0 curve on the parameter plane, it was ex- 
pected that the roots would move into the right half of the s-plane on the 
root locus, indicating instability of the loop. It was further expected 
that when the constant-alpha line crossed the £> = 1.0 line on the para- 
meter plane, that the real roots on the root locus would come together 
on the real axis and reach their breakaway point. The resulting root- 
locus plot is shown in Fig. 21 and it was found to fully verify the re- 
sults expected from parameter plane predictions. Additional roots which 
do not appear on this plot due to scale limitations were at -12.873+J3 .8547 
and -9.3502+jll.811. System zeros were at 0.0, 12.872+j3.8544 and 9.3502+ 
jll.811. As expected, the real roots which appeared at -1.2176 and -0.039322 
when beta had the value 0.1, moved together as beta increased. The left 
hand root was at -0.52 at (3 = 0.147 as predicted by the parameter plane 
and both real roots came together at 0.192 (^ = 1.0). The roots 

then became complex and moved toward the right-half plane, crossing over 
at (3 = 0.333 as predicted by the parameter plane. 

It was expected that closed-loop response could be improved by in- 
creasing beta, as predicted by the increase in damping ratio subsequently 
provided. The original system response to a step input is shown in Fig. 22 
with (3 55 0.147. Several other values of beta were tested and these re- 
sults appear in Fig. 23. To verify both root-locus and parameter plane 
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predictions of instability, a test with (3 = 0.333 was run and the re- 
sulting response is shown in Fig. 24. An almost undamped system was 
evident, characteristic of a system on the verge of instability with 
imaginary roots. 

On the basis of transient response tests, a new operating point of 
0< = 2.37, P-o .21 was chosen. Since the thermostatic level indicator 
transfer function information was somewhat doubtful as indicated in Ref. 

2, it was felt that this parameter could not be varied with any meaningful 
results. The newly selected operating point is marked with a G on the 
parameter plane. Fig. 20. 
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IV. SYSTEM RESPONSE TO VARIATIONS 



Using the complete simulation program once again, the system was 
subjected to the input shown in Fig. 5. This time, the parameters used 
on the parameter plane graphs of the previous section were varied to 
observe their effects on the overall system. Note that while a total 
of six parameters were investigated, it was recommended that only three 
of them be changed for "optimum" response. In the air and oil loops, 
the proportional gain of a controller was varied while in the water loop, 
one of the combining relay gains was changed. For ease of identification, 
these variables were subsequently called AIR, OIL, and WAT, respectively. 

It was decided to vary one parameter at a time first, to observe the 
effects, if any, of variable interaction. After this was accomplished, 
all three were changed in accordance with recommendations made in the 
previous section and the responses were again plotted. 

A. VARIATIONS IN AIR LOOP 

Table One indicates the values of the variable parameter AIR for 
the response curves of Figs. 25 through 29. 

TABLE ONE 



AIR CURVE # 

1.0 1 

1.67 2* 

2.0 3 

2.5 4 

3.0 5 

3.5 6 



^Original Operating Point 
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Curve # 4 was the result of using the operating point recommended 
from parameter plane predictions. It was a good deal faster than the 
original system as evidenced on the oil and air flow responses, Figs. 26 
and 27. It made little difference in the two most important variables, 
steam pressure and water level, Figs. 28 and 29. Note the near absence 
of coupling of the variation of AIR to water flow as evidenced in Fig. 25, 
the water flow response. The curves nearly coincided and were indistinguish- 
able as to which curve was which. 

B* VARIATIONS IN OIL LOOP 

Table Two indicates the value of the variable parameter OIL, for 
the response curves of Fig. 30 through 34. 



TABLE TWO 



OIL 


CURVE # 


4.0 


1 


5.0 


. 2 


6.67 


3* 


*Original Operating Point 





As was evidenced from the five response curves, Figs. 30 through 34, 
any proportional gain (OIL) less than 6.67 slowed the system responses 
considerably. This indicated that an increased gain, as predicted by 
the parameter plane analysis, was in order. 

C. VARIATIONS IN WATER LOOP 

Table Three shows the values of the parameter WAT for the curves 
depicted in Figs. 35 through 39. 
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TABLE THREE 



WAT 


CURVE # 


0.21 


1 


0.147 


2* 


0.1 


3 


0.05 


4 


★Original Operating Point 





All curves except those in Fig, 35 practically coincided with one 
another. Figure 35, the water flow rate response showed that the water 
responded faster with the new variable value (WAT = 0.21), as predicted 
on the parameter plane. 

D. VARIATIONS IN ALL THREE LOOPS 

As a final system test, all three parameters, AIR, OIL, and WAT 
were changed as was recommended in the previous section in accordance 
with parameter plane predictions. As a comparison, the old (original) 
response was plotted, along with a third curve representing the response 
if the parameters were moved in the opposite direction from that pre- 
dicted by the parameter plane. Table Four gives parameter values along 
with corresponding curve numbers for Figs. 40 through 44. 



TABLE FOUR 



AIR 


OIL 


WAT 


CURVE # 


2.5 


8.0 


0.21 . . . 




1.67 


6.67 


0.147 . . 


... 2 * 


1.0 


4.0 


0.1 .. . 





^Original Operating Point 
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In each case, it can be seen that the response is faster for curve 
#1 than for the original system. This proves the merits of design by the 
parameter plt-ne. The curve #3 in each case was considerably slower. 

Eo ANOTHER PARAMETER VARIATION 

It had been thought that the time delay (0.45 sec) in the feed water 
control valve contributed significantly to the oscillatory behavior of the 
water flow loop. The system was run using a variable time delay in the 
simulation to observe its effects. The time delays used were 0.45 (original) 
sec., 0.22 sec. (h ) , and 0.11 sec. . The response of the water flow 
loop is shown in Fig. 45. 

It can be seen that the variation in time delay did not affect the 
oscillations but did cause them to die out more quickly. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



From this work, it can be seen that the use of a simulation program, 
particularly DSL, was valuable in the testing and analysis of a system of 
this type. The value of other computerized analysis tools was also demon- 
strated. The parameter plane proved to be most valuable and, even where 
interaction was ignored, proved useful in the design of boiler control 
systems . 

The following recommendations are made as a result of this investi- 
gation: 

1. Apply the simulation technique to other boiler control 
systems, on other ships which have the common pneumatic variety and on 
the newer ships to come which will have a state-of-the-art (probably 
fluidic) system. 

2. Construct a three-dimensional parameter-space plot to 
depict the effects of varying the three parameters, OIL, AIR and WAT 
on the important variables, steam pressure and water level. 

3. Investigate the effects of varying other system para- 
meters or replacing certain system components, using the original simu- 
lation. 

4. Extend the original system to be dynamically sound over 

the entire operating range of the boiler when other transient response data 
becomes available. 

5. Continue work on the non-interacting controller shown in 
Appendix B. 

It is not the belief of the author that the improvements were of 
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sufficient magnitude to recommend changes to the actual system. It is 
hoped, however, that the methods and techniques used in this thesis will 
be of sufficient value to aid in the further design and improvement of 
boiler automatic combustion control systems. 
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APPENDIX A 



BOILER DESIGN AND PERFORMANCE DATA 

DESIGN 

Steam Generating Surface Area: Sq. Ft. 

Main tube bank 4975 

Screen tubes 933 

Front tubes . . 40 

Rear wall tubes 158 

Side and roof wall tubes 274 

Economizer 7500 

TOTAL WATER SURFACE 13880 

Superheater 2390 

TOTAL SURFACE AREA 16720 

Furance Volume: 1270 Cu. Ft. 

Boiler Weight: Lbs. 

Dry 230085 

Water (Steaming Level) . . . . „ 20080 

TOTAL BOILER WEIGHT 250165 
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PERFORMANCE DATA 



Operating Rate, 1 0 

Total Steam Generated, Ib/hr 

Drum Pressure, psig 

Superheater Outlet Pressure, psig 
Temperature, deg, F. 

Air Flow, Ib/hr 

Oil Flow, Ib/hr 

Anticipated Efficiency, % 

Number of Burners in Operation 



CRUISING 


FULL 


OVERLOAD 


44.7 


100 


120 


117478 


261450 


313800 


1210 


1245 


1265 


1200 


1200 


1200 


950 


970 


965 


145000 


321000 


386000 


8645 


20050 


24100 


87.6 


85.2 


84.2 


4 


6 


6 
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APPENDIX B 



A NON-INTERACTING BOILER CONTROLLER 

The utility of a proportional plus reset type controller is 
limited in an environment which has many variables in which inter- 
action is present. In the CVA-63 boiler control system, for example, 
a change in steam flow affects the drum water level as well as the 
pressure. When correction is applied to the water flow rate to bring 
the water level back to normal, the pressure is also affected which 
causes an adjustment to be made to the fuel flow rate. 

It is desired to design a controller which will automatically limit 
the interaction between control system variables. 

The major input variable to the boiler and control system is the 
steam flow rate. 

(la) AP d = Kps AWs + Kpf AWf + Kpi. AWc 

(ib) AX - Kxs AWs + Kx-fAW-f + Kxi AWl 

These equations determine the change in drum pressure ( APd ) and water 
level (AX ) in terms of steam flow rate change (AW s ) s fuel flow rate 
change ( AWf) y x nd feed water flow rate change (AWt). The coeffi- 
cients of Eqns. (la) and (lb) are time dependent and are determined from 
boiler dynamic analysis. 

The controller must adjust the fuel and water flow rates in a manner 
such th-.t when the st^am flow rate changes, neither the drum water level 
or the pressure is changed. 
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The controller will receive measurements of steam flow rate (W s ), 
drum pressure ( Fd ) and drum water level (X) from the boiler, compare 

/\ a 

the latter two with pre -determined set points, Fci and X and adjust flow 
of water and fuel accordingly. Refer to Fig. 46. 

Hence, the controller equations may be written: 

(2a) AWf ~ H-fp(Sp) + H-Fx (Gx) 4 Hfs (AWs) 

(2b) aWl - Hl p ( e P ) + Hu (ex) + Hu(aws) 

A A 

where <^p =p d-Pd and e*- X-X are error signals. The time dependent co- 
efficients Kcj are known (Fig. 46) as they represent boiler dynamic be- 
havior* The Hi coefficients are controller coefficients and are yet to 
be determined. 

If the equations (2a) and (2b) are substituted into (la) and (lb), 
the result is: 

(3a) A Pd ~ Kps AWs 4 Kpf [ Hfp (e p ) + H*fx + Hf5 AWs 1 

+ K p l [.Hip(ep) + Hix (e x ) *■ Hi* AWs] 



(3b) AX = Kxs AWs + Kxf fHfp (e P ) + H-f* (e*) + Hfs AWs] 

+ KxtlHip(6p) +• Hi.x (0x) + W is AWs ] 

Rearranging the above equations a.:d collecting terras, the above equa- 
tions can be written: 

(4 fl ) Z\P <3 = (Kps + Kpf Hf S + K P lHcs)aWs + (K P f FUp + Kpi Hep) e P 

(Kpf Ufx + Kpi 
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(4b) = (Kxs + K \4 \r\is 4 - Kxc WlS ) AWs + (Kxf H-fp 4 Kxc Hip ) e P 

+ (Kvcf W-Ck 4 Kxc H cx ) ex 

Requirements for the non-interacting controller force Af^ and AX 
to be zero for any change in €? or 6* # Therefore the coefficients 
of these terms must be zero identically. 



(5a) 


K P f Hf P 


4 


K pc Hep = O 


(5b) 


Kpf H-fx 


4 


K Hi-x ~ 0 


(5c) 


Kxf Hfp 


4 


Kxc Hip = O 


(5d) 


H-fx 


-f 


Kxc H Lx — 0 



Also present is the requirement t at when the pressure set point is changed, 
the pressure change must respond in some predetermined manner, say G(s). 

Then 

(5e) Kps 4 K p-f Hfs 4 Kpi. Hts - G p (S) 

Similarly, the water level must respond for a change in water level set 
point; therefore 

(5f) Kxs 4 KjC-F 4 Kxc Hi'S ^ GxCS) 

The six equations, (5a) through (5 f ) , contain six controller trans- 
fer functions ( Hej ) which are unknown, six boiler dynamic transfer 
functions which are known, and two n response to set point change" 
transfer functions which must be specified. Solving for the six un- 
knowns yields 
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(6a) i l-f s - (K* s Kpl - Kps Kxl )/ C 

(6b) Mts “ Kps - Kpf Kxs)/ C 



(6c) 


H f x 


- - (G* Kpu ) / C 




where 


C = K P f Kxi - K P 


(6 d) 


Hex 


= (G* KpO/C 


(6e) 


Hfp 


~ (.Gp Kxi ) /c 


(6f) 


H ip 


- -(GpKxO/C 



The simulation of this controller was attempted using DSL in the 
same fashion as the original control system. Several of the transfer 
functions contained more zeros than poles and the simulation failed. 

Future work on this controller could involve the introduction of 
excess poles in the offending transfer functions which could be cal- 
culated to minimize any disturbing effects on the system. 
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APPENDIX C 



I. PADE APPROXIMATIONS 

The general form of the fourth-order Pade approximation used in 
the simulations was: 

2 (rs) 4 - 20 (rs ) 3 4 i8Q(rs) 2 -84o rs + igso 

(rs) 4 + 2 o(rs) 3 + i8o(rs) 2 +840 rs + i68o 



Specific approximations were as follows for the individual loops: 



AIR: 



— . OST S 



e 






1.0551 x lo -5 S 4 - .003703 S + .58482 S - 47.8&S + 1880. 
(.08*51 *(0' s s 4 +■ , 00 3>7o'i> s 3 + ,58482 S 2 + 47.88 S + (680. 



OIL: 

- 025S . 39082 xio‘ fe 5 4 -. 3125 XI0~ 3 S 3 4- , 1125 S 2, - 2 1. S + 14,80. 

e • £ 

.39002 X (o'^s 4 4 ,3125 X|o' 3 S^ + .1125 s 2 + 21. S +1080. 



WATER: 



_.4SS 

e 



rv 



.04100625 S 4 - I.8Z25 S 5 ■(- 36-45 S 2 - 378. S + 1680. 
, 04 ( 004.25 s 4 + (.8225 S 3 + 36.45 S 2 + 376.S + K» 80 . 



II. INDIVIDUAL LOOP CHARACTERISTIC POLYNOMIALS: 

The following polynomials were calculated using the FORMAC com- 
puter program and were used in the parameter plane plotting routine: 
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AIR: 



.1240 x 10 " S ,s + .4B3I X lo -9 5" + .3625X I0 -7 s' 3 + . 8598 X lo“ 5 s‘ a 
+ .4735 x |o _i s" 4.01329 s'° 4 . 26135 S 9 4 3.5\34 S 8 ( 35.736 4 
.528?Xlo' 7 o<) S 7 + (258.293 - . 1691 x lo’ V +.52897 x lO~ 7 p) S% 

(l433. G + .002364 a - JG9IXIO -4 S 5 + ( 5264.03 15295CX 4. 002364 $ 4 
(14502. 64-1- 1.637 £3 - .15295(3 ) S 3 4 ( |2894. 84 + 204.6o3°(. + 1.637(3) S 2 
( I680. + 2056.32* + 204.603 p) S + 2056>,32(3 

OIL: .2973 x lo' T s ' 3 +.2611 X I0~ 4 s' 2 + .0|o3 s" + 2.144 s'° 

+ 197.64 S 9 + 1765.72 S 8 + 15266,9 S ? + ( 85501.7 + ,1104 X lo" & o( ) S 6 
4 ( 233383.7 - .96 X lo + .Uo4 X \o 5 ^jS + (286758.5+ .3756* - 96xlo~ 3 p ) S 4 
+ ( 123722. - 76.22 c< 4 , 3756(3) 5 3 4 ( |6SO. 4 66Z7.8 a “76.22 13) 5> 2 
4(29. 856 cX + 6627,8(3) S 4- 29.856 (S 

WATER: .03346 S 8 4 1.5422 S 7 4 32.24 S fo 4 (359.64 +, 127 XI0~ 4 (3 ) 5^ 

4(1924.32 4 ,6381 x |o - "cx -.5469 x |o" 3 ^') 5 4 4 (2748. - .283 x |o -9 ci 4 
.01129(9) S 3 +(2400.72 4 . 5672 x |o' 8 o( 4 .11718 (3 )s 2 + ( 1680. - 

.58824 x io" 7 * +.5208(3)5 4 . 2614 x IO -< ^0( 
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TITLE ***** CVA-63 BOILER AND ACC SYSTEM SIMULATION ******* 



***** INPUT FORCING FUNCTION ******************************* - 
AFGEN STEAM=0.0,C„ 0,6 3 0,0oG,8o0, 50000* ,ll n 0,0o0,100<> ,0 o 0 "" 

GS=AFGEN(STEAM,T IMF > 

********* **** + ** ********* **********<«********* *************** 
********** ********* ***• ******<t**tc^ *************>; *■ **** 



***** BOILER PORTION OF SIMULATION ************************* 

PDS= INTGRL ( ICl , - 1 9„ E-06*GS ) 

* PDF = TRNF R ( 1 . 2 , I C A , NUMA , DEN A , GF ) 

-POW=TRNFR( 0,3 ,'ICB, NUMB, DEN B , GW ) 

PQ=PDS + PDF +P0W 
PSA=1.37E-04*GS 
POUT =PO-PS A T 

LF = LEDL AG( IC2, 14o2, lo 4 2 , 5 a 1 1 E- 06*GF I 

LF1= INTGRL IlICULF) T T 7 

L W= I NTGRL (JC3,,5o02F-07*GW) " 

L SA=CMP X PL ( IC4, ; I C5, 2, 5 5,0*63, Oo 79E-05*GS ) 

LSB= INTGRL ( IC6 » 5» 0 2 E-0 7*GS )' 

LST=LSA-LSB 
L=LF1+LW + LST ^ 

*** £ $ j$r jfc * £ 4 9 ,:$ * * £ 4*;* 

# a? £ j*c a^c sjc j6c 4c 



***** OIL FLOW LOOP PORTION OF SIMULATION ****************** 
'PP=CMPXPLUC7,ICRj„0422f7o09,P0UT*1.25)' 1 ’ 

PZP=PRP-PPLV 

. PMA=PZP*0 1 L y / 

PMB= I NT GRL ( IC9 , 0 » l A*PZ P ) 7 
+PM=PMA+PMB / . tM 

S IGNAL=PQRL-PML 7 y 

PSF=FCNSW( S IGNAL,PQRL,PQRL, PML ) ' /-» 

PRF=CMPXPL( IC17, IC18,1«13, lo48,8oO*PSF) / 

PFINT=CMPXPL( ICl 9, IC20, 1 • 0 1 , 2« A85 , 36* 1*P RF L) x / 
vPF=DELAY(10,0»23,PFINT ) y 
GF = 117.*PF / ■ 



***** AIR FLOW LOOP PORTION OF SIMULATION ****************** 
PEQ=PML-PQRL^ 

POAA=A I R*P EG . 

POAB= INTGRL! I CIO , 0. 024* PEQ ) " 

POA=POAA+POAB ."■ -, / 

PQA=CMPXPL (^CJLLi4CI2j,0,5,7,91 ,QA*6o ) " 

PQR=TRNFR( 1 ,2 , I CC, NO! MC , DENC , POAL ) 

PQI = TRNFR( 1 ,2, fcr, NUNC, DENC, POAL ^ 

GB1 1 NT = C MPXPL ( ICl 3, ICLA,0. 2237 , 5.59, PQ1*53000. ) 

X GB1 = DELAY( 1 0, CC 05 i G B 1 1 N T I / 

Q81=CMPXPL ( ICl 5, I C 16, lo 534 , 0 o 42 1 , Oo 444E-03*G B 1 ) ^ 
QA=3„0*QB1^ 

****************************** ******************** ********** 
***, *************** ********************* ******************** 



***** WATER FLOW LOOP PORTLON OF SIMULATION 
/PSW=®22*PS A-WAT+PW A 
PDW=PRL+PSL-PSW < 

GWINT=TRNFR (0 ,3, I CJQJ* iNUJiQl* DEND,PDW)^" 
x GW = DELAY(10,0*45,GW ITTT ) . 



**************** 



PWA=0. 925E-C4*GWx 

PL=2.37*L 

PSL = PS + PL 

* + **>;;********* + ****** + •***** ***+****"*■********** + ***** ******** 
******* $*********. ********** *.*****#****** *********4 . ********** 
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** LIMITING FDR CERTAIN PNFU^A TIC SIGNALS ** 

*. PPL=LIMIT( -20, ,?0. , PP ) C/ 

PML=LIMIT( -20o » 2 0« » PM ) ' , 

P0RL=LIMIT (-2 0® »20o T PQR) 

PCAL = LIMIT (-2.0.C t20. , POA) 

PRFL = LI MIT(-60o0,60o t PRF » 

PQAL=LIMIT(-20 3 t20^,PQA) ^ 

************* ******** ****************** *•**•* ******** A********* 
A A- 4: ******* ****** *********** * ***. * * ** 4* V * * * * * * * <5 4 ** * * * *** ** 



***** CONTROL PORTION OF SIMULATION ************************ 
STOP AG ICAI2 ) t ICB<3 I , I CC ( 2) , ICD( 3) ,NUMA( 2 ) » DENA ( 3 ) ,NUMB< I ) , - 
DENB ( A ) t NUMC ( 2 ) ? DE NC ( 3 ) » NIJMD ( 1 ) » DEND ( A ) 

TARLE XCAJ 1-2 ) = 2^0® 0, ICB( l_-3 ) = 3*0^0, ICCJJ.=ZIs2*0*0 JL o«» * 

LC QXX-.X1 =3*0 oO,».; 

N UM A ( 1 - 2 ) =TC 6~2^r3o 1 8E-0A t DEN A ( 1 -3 ) =7 1 0 A , 1 o 0 , C . 0 , « • o 
NUMBI1 l=-A„ A2E-06,'DF.NB( 1-A) = 10 o 1,2o5AtK0t0o0,»o. 

NUMC < 1 -2 ) = <, 06A , 1 « 0 , DENC ( l- 3 J = 0 006A , „ 08 , 1 » 0 T o o • 
NUM0(l)=-3. lE0A,tDEND(l-A)=«816fl«5A32»ln 20A, UO 

C THESE PARAMETERS ARE FOR THE ORIGINAL SYSTEM 

PARAM 0IL=6o67 t A I R= 1 „ 67 , WAT = <, IA7 

PARAM PRP=0o0tPPL=0®0» PS = OoO 

CONTRL FI NT I M= 5 0* »DELT= o 05 

INTEG RKSFX 

END 

STOP 
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